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The properties of symmetric nuclear matter are investigated within the Green's functions ap- 
proach. We have implemented an iterative procedure allowing for a self-consistent evaluation of 
the single-particle and two-particle propagators. The in-medium scattering equation is solved for a 
realistic (non-separable) nucleon-nucleon interaction including both particle-particle and hole-hole 
propagation. The corresponding two-particle propagator is constructed explicitely from the single- 
particle spectral functions. Results are obtained for finite temperatures and an extrapolation to 
T=0 is presented. 

I. INTRODUCTION 

The evaluation of the saturation properties of nuclear matter from realistic models of the nucleon-nucleon (NN) 
interaction is one of the challenging testing grounds for many-body theories of quantum systemsQ, The strong 
short range and tensor components, which are required in realistic NN interactions to fit the NN scattering data lead 
to corresponding correlations in the nuclear wave function. The importance of these correlations is indicated by the 
observation that a simple Hartree-Fock or mean field calculation for nuclear matter at the empirical saturation density 
using such realistic NN interactions typically yields positive energies rather than the empirical value of -16 MeV per 
nucleon Q- 

While this argument on the importance of correlation effects in the nuclear wave function is based on a theoretical 
calculation only, more empirical evidence on these short range correlations can be deduced from the analysis of nucleon 
knock-out reactions0,Q. A recent analysis of the (e, e'p) reaction on 208 Pb covering a wide range of missing energies 
indicates that the occupation numbers for the deeply bound proton states are depleted by the same amount of about 
15 to 20 percent 4]. This depletion of the deeply bound hole states can be identified with the corresponding depletion 
of hole states in nuclear matter with momenta well below the Fermi momentum 0, 0, Q. The spectroscopic factors 
for these deep lying hole states should be determined by the tensor and short-range correlations mentioned above. 
Long-range correlations, on the other hand, lead to an additional reduction of the spectroscopic factors for states 
close to the Fermi surface Since these long-range correlations are sensitive to the collective excitation modes of 

the system, they should be different in finite nuclei as compared to the infinite system of nuclear matter. 
Various tools have been employed to account for correlations in the nuclear many-body wave function. These include 
the traditional approach, the Brueckner hole- line exnansion|lOj. and variational approaches using correlated basis 
functions^lj, L2|. Attempts have also been made to employ the technique of a self-consistent evaluation of Green's 
functions [l3L Il4| to the solution of the nuclear many-body problem. This method offers various advantages: (i) 
The single-particle Green's function contains detailed information about the spectral function, i.e. the distribution 
of single-particle strength, to be observed in nucleon knock-out experiments, as a function of missing energy and 
momentum, (ii) The method can be extended to finite temperatures, a feature which is of interest for the study of 
the nuclear properties in astrophysical environments, (iii) The Brueckner Hartree Fock (BHF) approximation, the 
approximation to the hole-line expansion which is commonly used, can be considered as a specific approximation 
within this scheme. 

Attempts have been made to start from the BHF approximation and include the effects of the hole-hole scattering 
terms in a perturbative wav |l5l Il6l Il7 | . For a consistent treatment, however, one should treat the propagation of 
particle-particle and hole-hole states in the in-medium scattering equation on the same footing. This turned out to be 
a rather ambitious aim. Starting from a single-particle propagator, which is characterized for each momentum k by 
one pole at the quasi-particle energy e(k), only, the in-medium scattering reduces to the Galitskii-Feynman approach. 
Trying to solve this equation one is confronted with the so-called pairing instabilitv[5l|l8l EES. 
These pairing effects can be taken into account by means of the BCS aDDroach [2lll22ll23| . At the empirical saturation 
density of symmetric nuclear matter the solution of the gap equation in the 3 Si — 3 D\ partial wave leads to an energy 
gap of around f MeV. Another approach is to consider an evaluation of the generalized ladder diagrams with "dressed" 
single-particle propagators. This means that the single-particle Green's functions are not approximated by one pole 
term but one tries to account for the complete spectral distribution. 

Attempts have been made to represent the spectral distribution in terms of three discrete poles 24] or in terms of four 
Gaussians|2^, [2(|. Indeed it turns out that such an improved representation of the single-particle Green's function 
leads to stable solutions. This finding is supported by the investigations of Bozek and Czerski pTl l28l I29I ] employing 
separable interaction models. 
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It is interesting to note that the same instabilities also occur in studies of finite nuclei|3fj, leading to divergent 
contributions to the binding energy from the generalized ring diagrams. These contributions remain finite if the 
single-particle propagators are dressed in a self-consistent way. 

In the present paper we want to present a method in which the equations for the one-body and two-body Green's 
functions for nucleons in nuclear matter are solved in a self-consistent way, keeping track of the complete spectral 
distribution in the single-particle Green's function. It turns out that the consideration of finite temperature helps to 
stabilize the numerical representation of the spectral distribution. Therefore we first determine the solution for finite 
temperature and then extrapolate to the case of T = 0. 

After this introduction we outline the formalism of the evaluation of Green's functions for many-body systems at 
finite temperature in Section 2. The results obtained for nuclear matter using the charge dependent Bonn potential 
CDBONN0| are presented in Section 3, where we also sketch some of the numerical details. Section 4 contains a 
short summary and the conclusions. 



II. GREEN'S FUNCTIONS 



In the Green's functions approach, physical quantities are expanded in terms of single-particle propagators. In a 
grand-canonical formulation, the one-particle Green's function can be defined for both real and imaginary times t, 

tr{exp[-/3(g - M AQro(x^t( x T)]} 
13(X ' ^ X ' 1 ] = tr{exp[-/3(g- M iV)]} ' (1) 

where T is the time ordering operator. It acts on a product of Heisenberg field operators ip{ x t) = e ltH ip(x)e~ ltH in 
such a way that the operator with the largest time argument t (or it in the case that t is imaginary) is put to the 
left. A minus sign is included for each commutation. The trace is to be taken over all states of the system with all 
particle numbers. /3 is the inverse temperature and /x is the chemical potential of the system. Depending whether 
t > if or t < if, the one-particle Green's function can be expressed by the correlation functions <7 > (x, t; x', if) or 
g l< (x, t; x', t'), respectively, where the time ordering in eq. has been carried out explicitely. Due to the invariance 
of the trace under cyclic permutations, it can be shown that the one-particle Green's functions obeys the following 
quasi-periodicity condition 

<?(x, t = 0; x', f ) = -e^ 5 (x, t = -i/3; x', if). (2) 

A hierarchy of relations defines the equations of motion for the Green's functions. The equation of motion for 
the one-particle Green's function involves the two-body potential as well as the two-particle Green's function 
<7n(x, £; • • • ;x"',t'"). In general, the equation of motion for the A^-particle propagator will be coupled to the 
(N + l)-particle propagator, if the Hamiltonian contains a two-body interaction. A good approximation scheme 
for <?(x, i;x',t') must be based upon an appropriate truncation for the two-particle propagator. Introducing the self 
energy S(x, t; x', if), the equation of motion or Dyson equation for the imaginary time one-particle Green's function 
reads 



.d_ 

dt 2m 



g(x, t; x', if) - J dt" J dx"S(x, t; x", t")g{x", t"; x', if) = 5(x - x')S{t - t'). (3) 



If the full two-particle propagator is replaced by an antisymmetric product of one-particle propagators, E(x, i;x',i') 
is just a Hartree-Fock self-energy 'E HF , which is real. In general, the self-energy will contain an additional complex 
contribution S c . We will consider the T matrix approximation for gu that contains all particle-particle and hole- hole 
ladders. 

In a translationally and rotationally invariant system in space and time, the correlation functions depend only on 
the difference variables x — x' and t — if . The real and positive spectral function may be defined using the Fourier 
transforms of the correlation functions along the real time axis 

A(k,u;)=g > (k,u J )+g<(k,u J ). (4) 

Since, in the Fourier space, the condition J5J can be transformed to 

g>(k,u) = e^-^g < (k,u) i (5) 



FIG. 1: The graphical representation of the T matrix 



the correlation functions become 

g>(k,w) = [l-f(u)]A(k,«i) (6) 
g<{k,w) = f(u)A{k,u), (7) 

where f(u>) = [e^^ - ^ + is the Fermi function. The coefficients of the Fourier sum, that takes into account the 
quasi-periodicity of g(x,t;x' ,t') (cf. eq. |(2Jl), can be expressed using the spectral function A(k,oj) 

g{k,z v ) = / . (8) 

J-oo 27r z v - ^ 

z v = -^Yp + /i are the (fermion) Matsubara frequencies with odd integers v. g(k, z u ) can be continued analytically to 
all non-real z. Using the Plemelj formula, the spectral function can be written as A(k, uj) = —2 lmg(k, lu + if]), where 
g(k, uj + irj) corresponds to the retarded propagator. 

By expanding the complex contribution to the self-energy in terms of one-particle Green's functions it can be demon- 
strated that it inherits all analytic properties of g. It it thus possible to write 

7T J_ 00 Z-LU 

Eq. is a prescription to determine the Green's function from the self-energy. In frequency-momentum space, the 
Dyson equation is an algebraic equation from which one can derive the spectral function to be 

ah \ -2Im£(fc,cj + ir}) . . 

A(k,u) = T2 — . (10) 

[w - £ - Re £(fc, oj)f + [Im E(fc, w + i V )} 2 

Because of the uj dependence of the self-energy, the spectral function has not quite a Lorentzian shape. The on-shell 
value of the real and positive quantity T(k,uj) = — 2ImS(fc, uj + irj) is nevertheless interpreted as the spectral width. 
The next step is to obtain the self energy in terms of the thermodynamic T matrix. This renormalized interaction 
takes care of the correlations induced by the strong short-range and tensor components of the nuclear two-body force. 
Graphically, the T matrix is depicted in Fig.^ Note that the arrows indicate both forward and backward propagating 
nucleons. The analytic structure of the T matrix can be deduced from a product of two Green's functions, so that 
T > and T < obey a similar boundary condition as the correlation functions 

T>(Q) = e^ Q - 2 ^T<(n), (11) 

Like the Green's function, the T matrix can be written in a spectral representation 

Using the Plemelj formula to separate the real and the imaginary part, one obtains from eqs. 1)1 111 and (|12ll 

T>(n)-T<(Q) = [ e ^ n ~ 2 ^ - 1]T<(Q) 

= -2ImT(n + i?7). (13) 
It is now possible to express ImS(fc,w + irf) in terms of the retarded T matrix [l3l IT^j 

ImE(fc,w + i?/) = -- [E > (k,oj)+T, < (k,w)] 

= ~lj Wfj x |^[< kk 1 T> ^ + ^|kk05 < (^c,0 + <kk'|T<(c, + c,0|kk0 3 > (^c,')] 
i r r+°o j / 

= 2 J j^rJ x ^( kk 'i Im n^+^+"/)ikk')[/(^)+^+^)]^(fc'^')- (14) 



4 



In the last line of eq. (|14f> . eqs. I©, J7J) and (|13|) have been used. The T matrix elements are anti-symmetrized. 
&(f2) = [e^ n ~ 2 ^ — is the Bose distribution function, which appears because hole-hole scattering diagrams are 
treated on the same footing as the particle-particle ladders in the T matrix approach. The pole in the Bose function 
&(f2) at f2 = 2fj, is exactly canceled by a zero in the T matrix 0,H2 such that the integrand remains finite as long as 
the T matrix does not acquire a pole. Such a pole may occur below a critical temperature Tc, a phenomenon which 
is often referred to as pairing instability. 

A useful assumption, especially for low densities, would be to allow for forward propagation in the T matrix equation 
only, since the phase space for the holes is very small. Then, T K = 0, the Bose function in eq. i|14|) disappears and so do 
the complications due to the pole in the T matrix at low temperature. The equations describing this approximations 
can be cast into the form of the BHF equations for finite temperature if one makes further simplifying assumptions 
for the spectral function A(k,u). 

The determination of the full T matrix requires the the knowledge of the product of two one-particle Green's functions 
with equal (imaginary) time arguments, as can be seen in the graphical representation in Fig. ^ 

9n(kuk 2 ,t-t') = g(k u t-t')g(k 2 ,t-t'). (15) 

The one-particle Green's functions can be expressed as Fourier series. Multiplication by e lZT ^ t ~ t ' and integration 
over the time variable from to — i/3 yields 

gi l {k l ,k 2 ,Z T ) = — ^2 g{kx,z v )g{k 2 ,Z T - z v ). (16) 



v— — oo 



Z T is the sum of two fermion frequencies, Z T = + 2/i, with even integers r. The spectral representation of the 
Matsubara Green's functions ca n be inserted into eq. ijPol) and the Matsubara sum is converted into a contour 
integration as described in Ref. The result is 

9 uh,k 2 ,z T )=r^ r^A {ku .)A { k 2 w- m - fi :'\ w 

J-oo 27r J-oo 2?r Z t - W - W' 

where the relation f(Z T —u)) = l — f(u>) has been applied. Expression (|17|) can be continued analytically. Substituting 
lu' = fi' — u, the real and the imaginary part of the retarded propagator can be separated (f2 real) 



on , o • s T 3 [ + °° dtt' ImgUk u k 2 ,n' + ir)) , r , . . 

gUkuk 2 ,n + ir)) = / - VIU n ' ' ^+ilm 9l ° 1 (Ai,*2,n + U7 , 



(18) 



where 

r + oo 

m n9^{ki £*n O -]- hi ) = — tt I 

2tt" 

In a partial wave expansion, the T matrix can now be determined as a solution of a one-dimensional integral equation 



/■ + 0O J 

lmgUki,k 2 ,n + ir,) = -7T — A(k u u;)A(k 2 , - w)[l - f(w) - /(O - u)]. (19) 



(q\T l J l , ST (P i n + ir,)\ q ') = ( q \V l l ST \ q ') + Y^ / d/c' fc' 2 (gl^^l/c') ^(i-, ^ + i?7 , fc') <fc'|T/,f7 + ( J P, O + . (20) 

Note that gfr(k±,k 2 ,fl + irj) has to be expressed in terms of the total momentum P and the relative momentum q 
of the particle pair, which requires an averaging over the angle between P = ^(ki + k2) and q = ^(ki — k2). This 
approximation procedure leads to a decoupling of partial waves with different angular momentum J. The numerical 
solution of eq. I|2l)[l enables us to use any nuclear two-body potential given in momentum space. The summation of 
the partial waves, 

(kk'|lmT(n + ir 7 )|kk') = i- (2J+l)(2T+l)( (Z (k,k')|ImT i f T (P(k,k'),0 + i ? 7)|g(k,k')), (21) 

(JST)l 

yields the T matrix in the form that is needed in eq. (|14|) . 

Finally, the Hartree-Fock contribution has to be added to the real part of £ 

Z HF (k) = ^ £ (2J+1)(2T+1) / ^<s(k,kO|^f ST |g(k,kO>n(A/), (22) 
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where n(k) is the density distribution 

—f(w)A(k,w). (23) 

-oo ^7T 

Note that this is a generalized Hartree-Fock contribution, since the full one-particle spectral function is applied. 
Eqs. HU, OH, (EH, (EH, ©j {HU j 03) an d OH provide a closed system of equations that have to be solved 

self-consistently. 

Before the numerical procedure that was applied to solve this system is discussed in more detail in the next Section, 
we will outline two different approximations to the full approach. 

One can think of a simplified set of equations, where the non-trivial spectral functions are replaced by the quasi-particlc 
expression 

A(k,u) = 2n 5(u - e(k)) (24) 

both in the two-particle propagator, g° v given by eq. (|T7J) and also in eq. ifHIl . This delta type spectral function is 
peaked at the quasi-particle energy and introduces an energy momentum relation to the model. The quasi-particle 
energy spectrum e(k) is derived from the following on-shell condition 

k 2 

e(k) = — +ReS(fe,e(fe)). (25) 
2m 

We will refer to this scheme as quasi-particle scheme. The reduced system of equations can be found in Refs. |32l|^3. 
Such calculations have been performed e.g. by Aim et al. for simple separable potentials at finite temperatures 19]. 
They find a critical temperature Tc, below which the system undergoes a phase transition to a superfluid state. Using 
realistic potentials, however, we found it very difficult to achieve convergence in the quasi-particle scheme even for 
T > Tq, since for a wide range of low momenta, eq. (|25|l has no unique solution, which indicates the limitations of 
the quasi-particle picture. This problem was also found in Refs. |l8l l34| . 

The BHF equations at finite temperature were formulated by Bloch and De Dominicis [35| . In this approximation, 
that was sketched after eq. OH , spectral functions of the quasi-particle type are applied, too. With respect to 
the quasi-particle scheme, the BHF scheme represents a further simplification, because it does not include backward 
propagation in the intermediate states. Since it allows for stable solutions, the BHF results are used as a reference 
for the full T matrix results. 



III. NUMERICAL DETAILS AND RESULTS 



The self-consistent solution is obtained in an iterative scheme. As a starting point, we use a quasi-particle T matrix 
in the first iteration. In the subsequent iteration cycles, g®j is constructed with the non-trivial spectral functions 
given by eq. IjlUfl. We fix the inverse temperature (3 and the chemical potential fi and do not change fi during the 
iterative cycle, which implies that the density of the system will vary unless self-consistency is achieved. The integral 
equation l|2(J|l is solved by a matrix inversion procedure. However, a pole subtraction as described in Ref. [3a ] is 
unnecessary, because g^ has no longer a quasi-particle form. Nevertheless, the integrand has to be sampled with some 
care in the vicinity of the quasi-particle peaks, since both the imaginary part and the real part vary rapidly there. 
We use between 40 and 100 integration points for the uncoupled partial waves. 

We have applied the CDBONN potential [3l| in our calculations. All partial waves up to J = 2 are included in 
the sum in eq. (|21fl . Higher partial waves do not contribute significantly to the off-shell structure of S. In contrast, 
the Hartree-Fock self energy includes partial waves up to J = 9. In Fig. El the thermodynamic T matrix in the 
quasi-particle scheme, i.e., using a quasi-particle g^ in eq. (|20|l . is compared to the shape of the full T matrix for 
zero relative momentum of the nucleon pair. For pair energies < 2e(P), the T matrix has no imaginary part in the 
quasi-particle approach. This becomes obvious if one looks at the quasi-particle approximation to <7j\, 

9u qp{P,V + ir/, q) = — — t-ttttZ 7TTT\ ^ ■ ■ 26 ) 
w it - (e(fci) + e{k 2 )) 6 + 177 

An imaginary part in the quasi-particle T matrix can only be formed if expression has a pole. For a nucleon 
pair with zero relative momentum, the sum of the single-particle energies in the denominator yields 2e(P). This leads 
to a sharp structure in the imaginary part of the T matrix in the quasi-particle approximation for Q > 2e(P). This 
structure is completely smeared out in the full T matrix calculation. In this case ImT is not vanishing also for large 
negative energies values of fl. 
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Once the T matrix is obtained on the (P, CI, g)-mesh, a three-dimensional interpolation has to be applied in order 
to carry out the transformation to the integration variables, (|k'|, 8' and a/), of the energy-momentum integrals 
in eq. (| 14|) . After the evaluation of the real part of the self-energy with the principal value integral in eq. JJJJj we 
interpolate the smooth functions ReS(fc,o;) and ImE(fc,cj) rather than the spectral function A(k,u>) to calculate 
Imj]]. The careful evaluation of the integral (|TT))l is one of the crucial steps of the self-consistent procedure. We 
actually consider 

/+oo 
—A{ki , lu + e(ki))A(ha , CI - e(h) - w)[l - /(w + e(fo)) - /(Q - e(fci) - w)], (27) 
-oo 27T 

where f2 = f2 + e(fci) + £(^2)- is the on-shell energy, which is defined in the same way as the quasi-particle energy, 
eq. 125|) . Note that the peaks of the spectral functions in eq. I|27l) are located around u> — and u> = CI, independent of 
ki and k 2 . This simplifies the construction of the integration mesh. Additionally, in the subsequent angle-averaging 
procedure, we can take advantage of the fact that Im g®j(ki, ki, CI) is always peaked around CI = 0. 
A sum rule for the two-particle spectral function was given in Ref. |37| for zero temperature. For the present case, 
this sum rule can be generalized to 

1 r +oa dCl 

/ — Im k 2 ,Cl + it]) = 1 - n(fci) - n(k 2 ). (28) 

This relation was used to check the numerical accuracy that was achieved for Im after performing the integration 
in eq. Q19f) . Mesh spacings and integration limits were adjusted such that both sides of eq. I|28|) do not deviate by 
more than 1% for single particle momenta up to k — 3000 MeV. 

After 6-10 iteration cycles, a self-consistent spectral function A(k,ui) is obtained. We have performed calculations for 
a range of chemical potentials between fi = —24 MeV and /i = 5 MeV at T — 10 MeV and for a range of temperatures 
between T = 5 MeV and T = 20 MeV at /i = —15 MeV. We found that it is difficult to perform stable and converging 
calculations below T = 5 MeV in our scheme. The reason for these difficulties is that the number of interpolation 
and integration mesh points has to be increased strongly in order to obtain stable results. The spectral functions for 
lower temperatures exhibit structures which require a treatment with a larger number of meshpoints. 
To extract information for temperatures below T — 5 MeV, we have extrapolated the retarded self-energy to lower 
temperatures, using the results of five stable calculations at T — 5,7,10,15 and 20 MeV. With the extrapolated 
self-energy, we have calculated the spectral function at temperatures below T — 5 MeV. The reliability of this 
extrapolation to lower temperatures will be discussed below. 

For the range of densities and temperatures we considered, we have found no signals of a pairing instability in the 
full calculation. The signature of this pairing instability would be a pole in the T matrix at CI — 2/i and zero 
total momentum of the nucleon pair, P = |ig. Although the pole appears on the real CI axis only at the critical 
temperature, Tc, the formation of the pole structure can be observed already at higher temperatures as a precursor 
effect 19]. Inspecting Fig. 2 we can see this precursor structure in the quasi-particle T matrix for P = 0. This 
structure is significantly reduced in the full T matrix. The depletion of the Fermi sea due to short-range correlations, 
in addition to the temperature-induced reduction of occupation at the Fermi surface, weakens the pairing correlations 
and no indications for a transition to superfluidity are found for the explored range of temperatures. Of course, we 
cannot exclude such a transition for lower temperatures and/or lower densities. 

Characteristic differences between the quasi-particle approximation and the full T-matrix approach can also be seen 
in FigUl which shows a comparison of the spectral functions for nucleons with momentum k = in nuclear matter of a 
density p around 0.3 fm -3 at temperature T — 10 MeV. While the quasi-particle spectral function tends to zero below 
u) — 11 = 200 MeV, the self-consistent result shows a large tail at negative energies. This redistribution of strength is 
related to the tail that was found in the imaginary part of the full T matrix in Fig [21 In the full calculation one does 
not constrain e.g. the construction of the self-energy at energies below the Fermi energy to the admixture of two-hole 
one-particle configurations, as it is done in the quasi-particle approach, but the self-consistent evaluation allows for 
general n-hole (n-l)-particle configurations. The double- hump structure in the quasi-particle spectral function reflects 
the fact that, for low lying states, eq. I|25|l has no unique solution (cf. paragraph after eq. The self-consistent 

spectral functions do not show this feature any more. The double hump can be interpreted as an indication for a 
strong coupling of the one-hole configuration to two-hole one-particle excitations. Obviously this coupling is reduced 
in the self-consistent calculation. 

Spectral functions for nuclear matter at different temperatures are displayed in Fig. 21 considering a value for /i = 
— 15 MeV, which corresponds to a density of about 0.36 fm~ 3 . The left panel presents the spectral function for low- 
lying hole states with momentum k = OMeV. Note that the spectral function for T = 0, which has been obtained by 
extrapolating the self-energy from self-consistent results evaluated for T larger than 5 MeV, shows the correct behavior 
and vanishes for to equal to fj,, separating the hole and particle part of the spectral function. This separation is of 
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course smeared out at finite temperatures. The width of the spectral distribution is large for k — at all temperatures 
and almost no broadening effect due to the temperature can be observed. This can also bee seen from Table [H which 
lists results for this width. 

The spectral function for loosely bound hole states is displayed in the right panel of Fig. 0] assuming a momentum of 
k = 345 MeV, which is just below the Fermi momentum of the density under consideration. In this case the width of 
the spectral distribution is much smaller than for k = 0. Note the different energy scale in this part of the figure. With 
increasing T, the peak of the spectral function clearly broadens and it is shifted to higher energies. At T = 10 MeV, it 
is almost symmetric around u) = p,. This means that at this temperature one does not observe any dip in the spectral 
function at the energy uj = p. 

The binding energy per particle, E/A was calculated from the self-consistent spectral function using the Koltun's sum 
rule: 



where d — 4 is the spin-isospin degeneracy factor for symmetric nuclear matter. A corresponding integral can be used 
to determine the kinetic energy per nucleon, Ekin/A, which implies that we can also evaluate the potential energy 
Ep t/A. The density p of the system is given by 



Both energy and density converged to a high degree of accuracy of typically better than 0.1% by the end of the 
iteration. The temperature dependence of some energy observables is given in Table IJ All values at T = 5 MeV and 
above were obtained from eq. (|29|l using self-consistent spectral functions. The values at T = 4 MeV and T — 3 MeV 
were calculated from the self-consistent results at T larger and equal to 5 MeV, extrapolating the results for the 
self-energies. The self-energy is a rather smooth function, but of course the extrapolation introduces uncertainties. 
We estimate an error up to about 3% for the interpolated values. In order to test this interpolation, we have also 
evaluated the total energy per nucleon for those temperatures, for which a self-consistent result exists, by interpolation 
from the remaining temperatures. These energies are listed in the column of Table which is labeled E/A ex ^ 
At low T, we found an almost perfect quadratic dependence between the binding energy E/A and the temperature. 
From this dependence we extrapolate the total energy per nucleon at zero temperature to -17.9 MeV . The decrease 
of binding with increasing T is mainly, but not only, due to an increase of the kinetic energy. 

In the last two columns of Tabled we present the on-shell width T(k,e(k)) at k = and at a momentum close to 
the Fermi surface. While this width is almost independent on the temperature for the low-lying states (k = 0), it 
increases strongly with T for the states close to the Fermi surface with a momentum (k — k* , see eq.(J30J)- 
The density dependence of the spectral width is shown in Table ITT1 The density of the system and the Fermi momentum 
k* of a filled sphere corresponding to that density are also reported. While the on-shell width for k = increases 
with density, the width at k = k* decreases. As we have seen already in Table U the width for the low-lying states 
is mainly determined by the phase space of the occupied states, which does not depend on the temperature, but on 
the density. In contrary, at the Fermi surface, the width is determined by the temperature deformation of the Fermi 
surface. Since this deformation is more drastic at lower densities, the width there is larger. 

In Fig. [SI the binding energy per particle E/A (solid line with squares) and the chemical potential p (dashed line with 
squares) are plotted versus the density of the system at T = 10 MeV. For a comparison we also present corresponding 
results of BHF calculations (continuous choice for the single-particle spectrum) for T = 10 MeV (thick lines) and 
for zero temperature (thin lines). We find a repulsive effect in the full T matrix calculation compared to the BHF 
calculation at the same temperature. This repulsive effect increases with density, therefore the saturation density 
obtained in the full T-matrix calculation is smaller than the one obtained in the BHF approach. The resulting 
density, p sa t — 0.31 fm~ 3 , however, is still almost twice as large as the 'empirical value', pp = 0.16 fm~ 3 . This 
repulsive effect relative to the BHF calculations is consistent with the observations of Ref . [27| , both with respect to 
the density dependence but also with respect to the size of the repulsion. 

The so-called Hugenholtz-Van Hove theorem states that, at zero temperature, the chemical potential should be equal to 
the binding energy at the saturation point [3^ ■ The theorem is clearly violated by about 20 MeV in the BHF approach, 
while it should be fulfilled in the self-consistent T matrix approach. This was first shown to be the case by Bozek 
and Czerski for a self-consistent T matrix approach using separable potentials [2J]. Although a zero temperature 
treatment is not really feasible in our approach, we can conclude a reasonable thermodynamic consistency in our 
calculation from the fact that the intersection point of p(p) with the binding curve at T — 10 MeV is located slightly 
above the saturation density. While the binding energy will increase as T approaches zero, the chemical potential will 
become less attractive for a given density, such that the intersection point moves closer to the expected region. 




(29) 




(30) 
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IV. CONCLUSIONS 



In this paper, we have presented a method for a self-consistent evaluation of the one- and two-body Greens function 
of nuclear matter in the ladder approximation. In contrast to previous works, no parameterizations of the single- 
particle spectral functions were used and and the full structure of a realistic NN interaction (CDBONN) was taken 
into account. This requires evaluation of the in-medium NN T matrix, the self-energy and the spectral functions of 
the nucleons for a wide range of energies and momenta. The equations describing these quantities were solved in an 
iterative scheme, for a given temperature and chemical potential. Self-consistency could be established after several 
iterations so that calculation of observables from these Greens functions converged to a high degree of accuracy. 
This computational scheme works rather well for temperatures above T — 5 MeV. The spectral functions for lower 
temperatures exhibit sharp structures, which requires a large number of meshpoints for a reliable representation. This 
inhibits direct calculations at very low temperatures. Therefore we introduce an extrapolation procedure, which is 
based on the smooth behavior of the nucleon self-energy, to deduce results also for lower temperatures. 
For the range of densities and temperatures we considered, no signals of a pairing instability have been observed in 
the full calculation. The distribution of strength in the self-consistent spectral function shields the system against this 
instability, which is observed in the quasi-particle approximation. This su ppo rts the findings of investigations using 
separable NN interactions or parameterizations of the spectral distribution|26l 111} • 

Comparing the total energies per nucleon calculated with the self-consistent T-matrix approach with the corresponding 
results obtained in BHF calculations we observe a repulsive effect of around 5 MeV per nucleon. As this repulsion 
increases with density, the saturation density of the self-consistent Greens function calculation is shifted to lower 
densities. The resulting saturation density, however, is still too large as compared to the empirical value. 
The consistent treatment of hole-hole and particle-particle scattering terms in the full T-matrix calculation leads to 
results, which are consistent from the thermodynamic point of view. Therefore, the Hugenholtz-Van Hove theorem, 
which is violated in the BHF approach, is respected in the Self-consistent Greens function approach (see also H3). 
We would like to acknowledge financial support from the Europaische Graduiertenkolleg Tubingen - Basel (DFG - 
SNF). 
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T [MeV] 


E/A [MeV] 


£/A ext . [MeV] 


E km /A [MeV] 


E pot /A [MeV] T(0, 


e(0)) [MeV] F(k* 


e(fc*)) [MeV] 


20 


-7.05 


-6.3 


58.41 


-65.46 


59.2 


15.5 


15 


-11.14 


-11.3 


54.45 


-65.61 


58.2 


9.3 


10 


-14.53 


-14.7 


52.48 


-67.01 


58.4 


4.7 


7 


-16.32 


-16.2 


51.87 


-68.17 


60.5 


2.5 


5 


-17.0 


-17.1 


50.9 


-67.9 


62.6 


1.5 


4 




-17.4 


50.9 


-68.3 






3 




-17.6 


51.0 


-68.6 







TABLE I: Energy, E/A, kinetic Energy, Ekin/A, and potential energy per nucleon, E po t/A, calculated for symmetric nuclear 
matter for a chemical potential p = —15 MeV, which corresponds to a density of p — 0.36 fm~ 3 , at different temperatures 
T. The values at T = 4 MeV and T — 3 MeV were obtained from an extrapolation of the self-energies evaluated at higher 
temperatures. Also displayed are the values for the on-shell width of the self-energy of nucleons with momentum k = and 
k = k* = 340 MeV (see eq.®). 



p [MeV] p [fm~ 3 ] k* [MeV] r(0,e(0)) [MeV] r(fc*,e(fc*)) [MeV] 



-23 


0.254 


307 


48.9 


6.7 


-20 


0.306 


326 


53.3 


5.5 


-15 


0.364 


345 


58.4 


4.7 


-5 


0.44 


368 


63.0 


3.5 


5 


0.51 


385 


70.1 


3.7 



TABLE II: Density-dependence of the on-shell width for k = and the momentum k* (see ea. 1301 '). The density is reported 
for each chemical potential, too. All calculations were performed at T — 10 MeV. 
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FIG. 2: Comparison between the quasi-particle approach (left panel) and our full scheme (right panel) on the level of the 
T-matrix. The dependence upon the cen " " ~ " " pair 

momentum q is set to zero. The chemical ] 



quasi-particle 

full approach 




I i I i i i I i I 
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FIG. 3: Spectral functions for nucleons with momentum k = at a temperature T = 10 MeV. The spectral function of the 
self-consistent calculation (solid line) has been determined for p, = —15 MeV, which corresponds to a density of p = 0.36 fm -3 . 
The results for the quasi-particle approximation (dashed line) are taken from a calculation at p = 0.3 fm -3 . 
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FIG. 4: Spectral functions in nuclear matter as a function of the energy at different temperatures. The left panel shows results 
for k — 0, while the spectral functions displayed in the right panel have been evaluated for k = 345 MeV, a momentum just 
below the Fermi momentum for the density p — 0.36 fm~ 3 , which corresponds to /i = —15 MeV. 
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FIG. 5: Binding energy per nucleon (solid lines) and the chemical potential (dashed lines) in symmetric nuclear matter as a 
function of the density p. Results of a full Green's function calculation at a temperature T = 10 MeV (squares connected by a 
line), are compared to BHF results are shown for zero temperature (thin lines) and T = 10 MeV (thick lines). The large square 
represents the 'empirical' saturation region that can be derived from the Bethe-Weizsacker formula. 



